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Abstract — This paper mainly modifies and further develops the comprehensive national strength model. By modifying the 
basic comprehensive national strength model, it can more accurately illustrate the society phenomena with time delay. First, 
we research the dynamics of the modified with time delay. By employing the normal form theory and center manifold method, 
we obtain some testable results on these issues. The conclusion confirms that a Hopf bifurcation occurs due to the existence 
of stability switches when the delay varies. Finally, some numerical simulations are given to illustrate the effectiveness of our 
result. 
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I. The Establishment of the Model 

Ordinary differential equation is used to study the comprehensive national strength of a country's overall situation, which is 
gradually developed in recent years. This article brings in time delay in the comprehensive national strength model. Besides, 
the comprehensive national strength model models are detailed analyzed, x = x{t) is hard power strength in / , which is a 
composite indicator the level of the material civilization (resource, economy, military, science and technology, etc.) of a 
country. As x(t ) bigger, the material civilization more prosperous. y(t) > 0 is soft power strength in t, which is the 
spiritual civilization of a country. y(t) > 0 stand for social evils (decision-making errors, education failure, official 
corruption people steal, etc.). At this time, soft power has obstacle to the social development. y(t) > 0 mean soft power is 
superior (wise decision-making, national quality, etc.), that has a promoting effect on social development. OC , [5 , y , m , M 
is constant. So, in the literature [1] comprehensive national strength shown by the following differential equation: 


,M-x s 0 
x = ax( ) - fly, 
M 

y = —yy + S(m - x)x. 


( 1 . 1 ) 


As is known to all, some change of the soft power can be reflected on the hard power after a certain period of time, the 
literature [2] has established the following delay model of the comprehensive national strength. 


* M — X 
x = ax( ) - J3y(t - r), 
M 

y = -yy + S(m- x)x. 


( 1 . 2 ) 


where is positive, and the other parameters are the same as of (1.1). While model of (1.2) had considered the time delay, in 
the real world, hard power strength not only rely on hard power strength of soft power in the past, but also rely on soft power. 
In order to make the model more accurate, we modify the system of (1.2) to the following form: 


• M — X 

x = ax( — ) - J3(ky(t - z) + dy), 

M 

y = - yy + 8{m - x)x. 


( 1 . 3 ) 
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where k + d = 1 . The organization of this paper is as follows: regarding T as bifurcation parameter, we study the stability 
of the equilibrium point of the system (1.3) and Hopf bifurcation of the equilibrium depending on r . Then, based on the new 
normal form of the differential-algebraic system introduced by Chen et al. [3] and the normal form approach theory and 
center manifold theory introduced by Hassard et al. [4], we derive the formula for determining the properties of Hopf 
bifurcation of the system in the third section. Numerical simulations aimed at justifying the theoretical analysis will be 
reported in Section 4. Finally, this paper ends with a discussion. 

II. Stability and Local Hopf Bifurcation Analysis 

The stability and Hopf conclusions for the system of (1.1) can be obtained directly from the literature [5,6]. From system 
(1.3), we can see that there exists equilibriums. 


E 0 =(X q ,Y q ) = ( 


M ( ay — mj3S) Mad ( m — M ) ( ay — m J3S) 


ay - M(d8 


(, ay-M{38 ) 2 


),^=( 0 , 0 ). 


According to the practical significance of the model, here we only discuss the problems of the hopf bifurcation and stability 
for the sole positive equilibrium point E 0 . 


Set U x (t) = x(t ) — X 0 ,U 2 (t) = y(t)-Y 0 . The system (1.3) becomes 


u x (t) = a{ 1 + p)u x w, 2 - /3(ku n (t -z) + du*,), 

< M “ “ ( 2 . 1 ) 

u 2 (t ) = -yu 2 + S(m + Mp)u [ -Suf. 


Where p = 


2 (ay - m fid) 
-ay + MJ35 


.And, the linearization of system (2.1) at E 0 is 


Mj(0 = a( 1 + p)u x - J3(ku 2 (t -r) + du 2 ), 

i (2.2) 

u 2 (t) = -yu 2 + S(m + Mp)u x . 

The characteristic equation of system (2.2) is as follows. 

— c/ i — ge + f — 0, (2.3) 

where c = a(l + p) + S{tn + Mp), f = aS(l + p){m + Mp) — dy/3, g = y/3k. To study the stability of the 
equilibrium point E 0 and bifurcation system, we only need to discuss the distribution of the roots for characteristic Eq. (2.3). 
If the all roots of equation (2.3) have negative real part, the equilibrium point E 0 is steady. If the equation has a root that 
contains positive real part, the equilibrium point E 0 is not stable. In order to study the distribution of the roots for Eq. (2.3). 
Considering T = 0 in first, characteristic equation (2.3) is 

At — cA, + q = 0. (2.4) 

Where q = aS( 1 + p)(m + Mp) - yfi . 

(H l )c<0,q>0. 
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Obviously, all roots of equation of (2.3) has negative real part if (H l ) is satisfied. So the system in equilibrium E 0 is 
locally asymptotically stable for T = 0 . Now, we investigate the local stability around the positive equilibrium point for the 
system (1.3) and the existence of Hopf bifurcations occurring at the equilibrium point E () when T 0 > 0 . 


Lemma 2.1 For the system (1.3), Eq. (2.2) has a part of purely imaginary roots for T 0 > 0 when (77,) is satisfied. 
Proof. If X = ico is a solution of the characteristic equation (2.2), when and only when CO meet 

-or - ccoi - g{coscoz - i sin cot) + f = 0. 

The separation of the real and imaginary parts yields. 


j-co 2 + f = gCOSCOT, 
[cco = gsincoz, 


(2.5) 


which lead to 


v 2 +(c 2 - 2f)v + f 2 - g 2 =0, (2.6) 

2 

where V = C0~ . We assume that the coefficients satisfies the following conditions. 

(H 2l )c 2 -2f < 0,4(/ 2 - g 2 ) < (c 2 -If) 2 ,{H 22 )c 2 -2f > 0,4(/ 2 - g 2 ) <0. 

If the condition (//-, j) or ( H _, ) is satisfied, then Eq. (2.6) has positive roots. Therefore, Eq. (2.2) have purely imaginary 
roots. From (2.6) we obtain 


CO Q =- 


l2f-c 2 +y [(^2f) 2 -4(f 2 -g 2 ) 


T J = 


— (arccos (— ) + 2jn), j = 0,1,2,- • • . 

8 


As a result, when T = Tj , the characteristic equation have a pair of purely imaginary root. The proof is completed. 

Lemma 2.2 The transversality conditions Re( — — ) 1 > 0 or ( H,, ) is satisfied. 

d r n 


Proof. By differentiating both sides of Eq. (2.3) with regard to T and solving X(t) . We have 


dX _ - gXe Ar 

cIt 2X-c + gTe~ Ir 


then 


Re A" 1 C' =Re( 2/t 

d t J - gXe 1 
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- 2co q costy 0 r ; + c sin o 0 r . 
g ® b 


c 2 -2f + 2a> 2 
g 2 


d^i 

We can get C 2 — 2/ + 26U > 0 by condition (H 2l ) or (H^) , thus Re( ) 1 > 0 . The proof is completed. 


Lemma 2.3 For Eq. (2.3), if T e [0, f 0 ), all of his roots have negative real part. Positive equilibrium is asymptotically 
stable, and the positive equilibrium produces Hopf bifurcation in T = T () . 

III. The Direction of The Hopf Bifurcation and The Stability of Periodic Solutions 

In §2, we obtain the conditions of Hopf bifurcation. In this section, we discuss the direction of Hopf bifurcation and the 
stability of the bifurcating periodic solutions based on the normal form introduced by Chen et al. [3] and the center manifold 
theory introduced by Hassard et al. [4], 

In the following part, we assume that the system (1.3) undergoes Hopf bifurcation at the equilibrium E 0 for T = T . , and we 
let i C0 Q is the corresponding purely imaginary root of the characteristic equation at the equilibrium E 0 . Set 
T = T q + jU, /U e R , clearly, /J — 0 is the Hopf bifurcation value of system (3.2). Set t = ST , Ui(t) = Uj(tT),i = 1,2 , for 
convenience, we continue to use ll^t) said lli(t) .Then the system (1.3) is equivalent to the following Functional 
Differential Equation (FDE) system in C e C([— 1,0], R~) , 


u(t) = L +F(ju, u t ). 


(3.1) 


where u(t) = (w 1 (f),M 2 (f)) r , L C — > R, F : RxC — > R are given, respectively, by 


L u (<f>) = (T 0 + W(0) + (r 0 + ]U)B 2 <t>(- 1), F(jU, <f>) = (r 0 + ju) 


v -^2 2 (°) J 


where L (j) = {(j) l ,(/) 2 ) T e C([— 1,0], R 2 ) . 


*1 = 


' a( 1+ p ) 

-pd ' 

= 

i 

o 

v -r 

8(m + Mp) j 

7 z 

10 0 j 


By the Riesz representation theorem, there exists a matrix function whose components are bounded variation function 
rj(6,ju) : [—1,0] — > R 2x2 such that 


LJ = \_d r/(0, ju)(f){0), f)eC. 


(3.2) 


In fact, we take 
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q{9, ju) = B { 8{9) + B 2 S(9 + 1). 

S{d) is a Delta function. For (j) e C ([—1,0], ' ) , we define operator A and R as follows. 

'dim) 


(3.3) 


A{R)m= 


d0 


^e[-l,0), 


\[di V id^)m\ 0 = 0 . 


(3.4) 


Roam = 


fo. 


0 e [—1,0), 


9 = 0. 

Then the system (2.1) can be written as the following form: 

u t = A(ju)u r + R(ju)u t 


(3.5) 


(3.6) 


Setting (p eC ([0,1], (R )*), the adjoint operator A of A is defined ; 


A*y/{s) = 


di//(s) 

ds 


s e (0,1], 


(3.7) 


[J ' _d(rf (t,0)y/{-t)),s = 0, 

and a bilinear inner product is given by 

(p{s),(/>{9) >= (p(0)(/>{0) - f° f (pi% - 9)drj{9)<j)i%)d p. 

J- 1 Js=0 


< 


(3.8) 


From the discussions in Section 2, we know that i TjQ) 0 are eigenvalues of A . Thus, they are also eigenvalues of A" . Next 
we calculate the eigenvector q(9) of A belonging to iTjCO 0 and eigenvector q* (s) of A * belonging to the eigenvalue 

-iT-(O 0 . By the definition of A , we have A(0)q(0) = — — — - and q(9) = q(Q)e' r, "' n ° [ n addition. 


d9 


J dq(9)q(9) = B x q{0) + B 2 q{-\) = A(0)g(0) = io) 0 q( 0), 


(3.9) 


Let q( 0) = (1, q 2 ) T , Eq.(3.9) becomes 


a(l + p) -pid + ke 0i ) 
v -y S(m + Mp ) 


ia>„Ti Af J \ 


\fl2j 


= ICO, 


f 1 \ 


Vh) 


where q 2 = 


a{ 1 + p)-ia> 0 
p(d + ke 


By the definition of A*, we have A*(0)<7*(s) = — and q* (s) = q* {0)e' TaCO<>s . Then, 

ds 
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| q* (-t)dq(t) = B^q*{ 0) + B T 2 q* (- 1) = A*q*( 0) = -ico 0 q*(0). (3.10) 

Let q*( 0) = D(l,q 2 ) T - Eq.(3.10) becomes 


a(l + p ) -y ' 

m 

. * 

rr 

K - J3(d + ke u "" Tj ) S(m + Mp) y 

* 

\ c h) 

= -ico Q 

* 

\di) 


* a(l + p) + ia> 0 
where q 2 = 

r 


<q\s),q(9)> = q (0)<?(0)-f f q (2; - 9)dr/(9)q(Z)d2;. 

J-l Js=0 

= D[l + q 2 ^ 2 * -\jl,Z)9e~ iT ^ 9 dq(9)( 1 1 ] 


= D[l + q 2 q 2 + t e " 0 ° r> (l,q 2 )B 


Vh) 

ro 

V' 


<h ) 


= D[\ + q 2 q 2 - J3kv jq 2 e " A>Tl ] , 
where (q* (s),q{9)) = 1 , soD = [1 + q 2 q 2 - fdkz fl 2 e ia> ° T >y l . 

Next, we study the specific parameters of the direction and size for bifurcated periodic solutions. Using the same notations as 
in Hassard et al. [4], we first compute the coordinates to describe the center manifold C (l at // = 0 . Define 

z(t) =< q*,u t >,W(t,9) = u t {9)-2Re{z(t)q(9)}. (3.11) 

On the center manifold C 0 ,we have 

W(t,9) = W(z(t),~z(t),9) (3.12) 

2 _ ; 2 

Where W(z, Z, 9) = W 20 (9) *-+W n (0)zZ + W m (0) y + • • • . 

In fact, z and z are local coordinates for center manifold C () in the direction of q and q . Note that W is real if U t is 
real. We consider only real solutions. For the solution of Eq.(3.6) in center manifold C Q , we have 

z(t) = (q*(s),jU t ) = (q * (6»),(A( 0) + 

= (A*q\s),ju t ) + q (s)R(0)ju, - [ f q -9)dq(9)A(0)R(0)ju,(^)d^ 

J-l Js = 0 


= (i z-co^i ( s), p t ) + q (0)F(0, p, ( 9 )) 


(3.13) 
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Let 


where 


z(t) = ico Q z + g(z,z). 


2 _ ~ 2 

g(z, z ) = g 20 — + g n zz + g 02 — + • 


from (3.13) and (3.15), we have 


— — * 


g(z,z) = q ( 0 )F 0 (z,z) = tV 


r n 


— * 

V^2y 



= T j V(F l +q 2 F 2 ). 


By (3.11), we have 


where 


u t (9) = (Mj (t + 9), u 2 (t + 9)) = W(t,9) + z(t)q(9) + z(t)q(0) 


M, (t + 0)= W (1) (t, 9) + z(t)q w ( 9 ) + z(t)q ' (9) , 


(i)/ 


- . I 1 ) 


u 2 (t + 9) = W (2 \t,9) + z(t)q u \9) + z(t)q (9) . 


( 2 ). 


- “( 2 ). 


Then 


g(z, z) = Tj - Sq 2 ql )z 2 + ^ - 28q 2 q 2 ql )zz + -8q 2 q 2 )z + 


M 


M 


M 


(2W 1 ( i 1) (0) + ^ (0) + ^ W 2 ( 0 2) (0)) - 8 q 2 (2W™ (0)q 2 + W 20 2) (0)<7 2 )]z 2 z + • 


Comparing the coefficients with (3.16), it follows that 

g 20 =2Tj D (-^--8q 2 q;), 
M 


#11 =' r J £> (-^- 2 ^ 2 ^ 2 2 ). 
M 


S m = 2z j F>{-^--8q 2 q 1 2 ), 
M 


g 2t = 2r . (2<> (0) + i (0) + ± (0)) - (2<> (0)$ 2 + W 2 ( 0 2) (0) ?2 )] . 


(3.14) 


(3.15) 


(3.16) 


(3.17) 


By Eq.(3.15), (3.16) and definition (3.1 1) of W(t, 9) , we have 
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W = u-zq-zq = 


A(0)W-2Rzq (0 )F o (z(t),z(t))q(0), 0e[-l,O), 

A(0)W - 2 Re ~q (0 )F 0 (z(t)~z(t))q(9) + F 0 , 9 = 0. 


=AW + H(z,z,9) 


Let W = A(0)W + H(z, Z.ff) . Where 


_ 2 _ - 2 
H(z, z,9) = H 20 (9)^- + H u (9)zz + H 02 (9)^- + - 


From the definition (3. 1 1) of W(t, 9) , we have 


W = VF_ z+W-z 


By (3.18), (3.19) and (3.20), we have 


H 20 (9) = 2i t 0 coW 20 (9)-A(0)W 2 j 0 (9), 
H n (9) = -A(0)W n (0). 


For 9 G [—1,0] , by (3.14), (3.18) and (3.19), we get 


H (z, z,0) = -2Re[q (0 )F 0 q(9)] 


-2 


= -(^20 y + 8 uzz + <?02 y + 821 y-)tf(0) 


_ 2 _ 2 _ 2 _ 

~(8 20 y + 8 uzz + g 02 y + <? 21 y-)tf(0) 


Comparing the coefficients with (3.19), we gives that 


H 20 (9) = -g 20 q(9)-g 02 q(9), 
H n (9) = -g n q(9)-g n q(9) . 


It follows from (3.21) that 


W 20 (9) = 2 icor 0 W 20 (9) - A(0)g 20 q(9 ) + g 20 q(9), 
Wu(9) = g u q(9) + g u q(9). 


(3.18) 


(3.19) 


(3.20) 


(3.21) 


(3.22) 


(3.23) 


(3.24) 


Solving for W~, {) (9) and W u ( 9) , we obtain 
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W 20 (6>) = l ^q(0)e iToO,e + ^^q(O)e- iToC00 
T a co 3t q co 

W x , (9) = - Mil q ( 0)e i(OT ° e + Mil q(0) e - ia,T ° e 
X 0 G) TqCO 


+ E x e 2iTo<oe , 

+ E 2 . 


(3.25) 


By (3.4) and (3.21), we get 

J IdqiOWje) = 2i Tj(o 0 W 20 (0) -H 20 (0) , 

\\dq{G)W n {G) = -H n { 0) 

From (3.17) and (3.18), we get 

(0) = ~g 2 (M°) - 8 02 9(0) + 2r y (£, , , £ 12 , £ 13 ) r , 


# i t(<9) = —^11^(0)— ^,i $(0) + r ; . (£ 21 , E 22 ,E 2i ) T . 

According to Eq.(3.25), (3.26) and (3.27), we get 

(2/r.® 0 /-J \e 2iT ^dqmE x =lr j (E u ,E n ,E u ) T , 


Then 


£, = 2 | 


2z<» 0 - a(l + p) P(d + ke ,m ° Ti ) 
y 2ioj () - S(m + Mp) 


y 1 

( -a ' 


M 

) 



e 2 =- I 


^a(l + p) “/? 

— y S(m + Mp) 


y 1 

( -2a ^ 


2 _ 

) 

2Sq 2 q 2 j 


According to the above Proposition and [7,1 1], we can compute the following parameters: 

C l (0) = -^(*,„* ll -2|g M |-i| ) ; 0; | 2 ) + 1 ? i , 

2z-j<» 0 3 2 

RejQCO)} 

2 ReR'^)}’ 

Pi = 2Re{C 1 (0)}, 

T _ Im{C 1 (0)} + // 2 (ImU'(r j )}) 

^2 — 

O) 0 Xj 

IV. Numerical Simulation 


(3.26) 


(3.27) 


(3.28) 


Numerical simulation shows from stable to unstable complex transformation process the system (1.3). We give a concrete 
example to show the dynamic behavior of comprehensive national strength model. We 

take OC — 0. \,M — 10, [3 = 0.3, k = 0.8, d — 0.2, y = 1 8 — 2,m = 9 in the system (1.3). So we can obtain 
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equilibrium E 0 =(11.8,12.93). Through calculating of the Mathematical software, we obtain C0 Q ~ 2.62 ,T 0 ~0.88. By 
lemma 2.2 and above conclusion, we can obtain that equilibrium E 0 is stable if T = 0.7 < T {) , (see Figure 1). By contrast, 
the equilibrium E 0 is unstable if r = l.l>r 0 , (see Figure 2). Besides, whenr 0 =0.886, the periodic solutions occur 
from the equilibrium E 0 (see Figure 3-4). However our analysis indicates that the dynamics of the comprehensive national 

strength model with time delay can be much more complicated than we may have expected. It is still interesting and inspiring 
to research. 



Figure 1 The equilibrium E 0 
OF SYSTEM (1.3) IS STABLE WITH T = 0.7 < Z 0 


Figure 2 The equilibrium E 0 
OF SYSTEM (1.3) IS STABLE WITH T = 1 . 1 > T 0 


y LdL 



Figure 3 System (1.3) produce the periodic solutions with t 0 = 0.886 
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Figure 4 Time sequence diagram of Figure 3 
V. Conclusion 

These papers apply delay in the comprehensive national strength model, which show rich dynamics behavior. Different from 
previous studies, we added the influence of time delay feedback in the system (1.3). Dynamic behavior of comprehensive 
national strength model with time delay is analyzed by using the method of quantitative. When delay $\tau$ across a series of 
critical value, nonlinear dynamic system generate the Hopf bifurcation. In addition, by employing the normal form theory 
and center manifold method, we obtain some testable results on these issues. We use normative theory and center of popular 
theorem obtained the calculating method of the direction of Hopf bifurcation and stability of periodic solutions. Finally, the 
above theoretical analysis is verified by numerical simulation. The dynamic behavior of the comprehensive national strength 
model is rich. Many aspects is not mining, yet to be. References 
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